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The energy positions of frontier orbitals in organic electronic materials are often studied experi¬ 
mentally by (inverse) photoemission spectroscopy and theoretically within density functional theory. 
However, standard exchange-correlation functionals often result in too small fundamental gaps, may 
lead to wrong orbital energy ordering, and do not capture polarization-induced gap renormaliza¬ 
tion. Here, we examine these issues and a strategy for overcoming them by studying the gas phase 
and bulk electronic structure of the organic molecule quinacridone (5Q), a promising material with 
many interesting properties for organic devices. Experimentally, we perform angle-resolved pho¬ 
toemission spectroscopy (ARUPS) on thin films of the crystalline /3-phase of 5Q. Theoretically, we 
employ an optimally-tuned range-separated hybrid functional (OT-RSH) within density functional 
theory. For the gas phase molecule, our OT-RSH result for the ionization potential (IP) represents 
a substantial improvement over the semi-local PBE and the PBEO hybrid functional results, pro¬ 
ducing an IP in quantitative agreement with experiment. For the bulk crystal, we take into account 
the correct screening in the bulk, using the recently developed OT-SRSH approach, while retaining 
the optimally-tuned parameters for the range-separation and the short-range Fock exchange. This 
leads to a band gap narrowing due to polarization effects and results in a valence band spectrum 
in excellent agreement with experimental ARUPS data, with respect to both peak positions and 
heights. Finally, full-frequency Go Wo results based on a hybrid functional starting point are shown 
to agree with the OT-SRSH approach, improving substantially on the PBE-starting point. 


I. INTRODUCTION 

Organic semiconducting devices are extensively stud¬ 
ied as they have the advantage, compared to their 
Si-based counterparts, of being flexible, cheap, light- 
weighted and processable at low temperatures.®® Re¬ 
cently, the organic pigment 5,12-dihydro-quino[2,3- 
b]acridine-7,14-dione (quinacridone), C 2 oH 12 N 2 02, a 
derivative of pentacene, has attracted considerable in¬ 
terest as an active material in organic electronics,®® e.g. 
due to its remarkable air-stability, high photogeneration 
efficiency, electrochemical stability, and high fluorescence 
lifetime in solution. 

Quinacridone (5Q) differs from pentacene by the pres¬ 
ence of functional N-H and C-0 groups at phenyl rings 
2 and 4 (see Fig. la). On the one hand, these polar 
groups cause the formation of intermolecular hydrogen 
bonds. These are responsible for 5Q’s ability to form 
self-assembled, supramolecular structures and allow for 
the remarkable air-stability of 5Q devices reported in the 
literature.® However, these functional groups also break 
the conjugation of the molecule and thereby significantly 
change the energetic positions and spatial shapes of the 
frontier orbitals. The latter are crucial for the perfor¬ 
mance of organic semiconducting devices because they 
determine electron removal and insertion energies and 
therefore the band gap and level alignment. As such, they 
are extensively studied in the context of organic molec¬ 
ular systems and organic/inorganic interfaces.®® 2 Such 
investigations strongly benefit from accurate and com¬ 


putationally inexpensive theoretical models, which are 
often necessary for a better interpretation of experimen¬ 
tal resultsP® 

An efficient work horse for first principles calcula¬ 
tions of the electr onic structure is density functional 
theory (D FTYp ^® usually within the Kohn-Sham (KS) 
framework!^ 61 ' In this scheme, the original many- 
electron problem is mapped uniquely into a fictitious 
noninteracting electron system, yielding the same elec¬ 
tron density. This mapping leads to effective single¬ 
particle equations that provide a significant conceptual 
and computational simplification of the original many- 
electron problem. However, due to the fictitious nature 
of the noninteracting electrons, the correspondence of KS 
eigenvalues with ionization energies measured in an ex¬ 
periment is not straightforward. 18 One exception is the 
highest occupied orbital whose energy can be rigorously 
identified with the ionization potential (IP) of the neutral 
system, if the exact exchange-correlation (xc) functional 
is usedPMUl j n g enera p Rg results based on approximate 
xc expressions, e.g. local or semi-local functionals, may 
suffer from pronounced self-interaction errors (SIE) and 
a lack of the derivative discontinuity in the xc poten¬ 
tial, and therefore are not expected to agree with ex¬ 
perimental findings obtained, e.g. from photoemission 
spectroscopy.®^ These conceptual problems may, at least 
partly, be cured by introducing the generalized Kohn- 
Sham (GKS) scheme,and considering hybrid function¬ 
als with a fraction of exact exchange either globally or in 
a range-separated manner. 23 ^ However, band gaps, IPs, 
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electron affinities (EAs), and the orbital order obtained 
within such hybrid functional may still be in error when 
compared to experiment! 25 * 26 * 

A promising strategy to improve the agreement with 
experiment is the more recent class of range-separated 
hybrid (RSH) DFT functionals,* 27 * 28 * where the interelec¬ 
tron Coulomb repulsion term is separated into long-range 
(LR) and short-range (SR) components. The LR term is 
mapped using Hartree-Fock theory, thereby establishing 
the correct asymptotic potential. The SR term is mapped 
using a (semi-)local KS functional or a conventional hy¬ 
brid functional, maintaining the compatibility between 
the exchange and correlation expressions. In this ap¬ 
proach, one still needs to determine the range-separation 
parameter as a universal value usually leads to en¬ 
ergy levels that, although greatly improve the accuracy 
of standard hybrids, are still not at the desired accu¬ 
racy level 12M30] This can be impro ved b y using optimally 
tuned RSH (OT-RSH) functionals,^ 5 31: where the range- 
separation parameter is tuned for each system such that 
physically motivated tuning conditions are fulfilled with¬ 
out introducing any empirical parameters. DFT calcula¬ 
tions using OT-RSH functionals have been shown to pro¬ 
vide an accurate, non-empirical description of band gaps, 
IPs and EAs for a variety of systems, among them atoms, 
molecules and polymer^2H251 as well as larger gas phase 
organic molecules relevant for organic semiconducting 
devices .®£ 3 Z j n addition, it has been demonstrated that 
the description of deeper lying occupied states of an iso¬ 
lated molecule can be considerably improved within the 
OT-RSH approach, if an additional degree of freedom 
is introduced, which allows the modification of exact 
exchange in the short rangePSEH In particular, it was 
shown to well-describe the outer-valence spectra of sev¬ 
eral organic molecules that exhibit a mixture of localized 
(tr) and delocalized ( 7 r) states, O* a challenging situation 
where the difference in self interaction error (SIE) for dif¬ 
ferent orbital types can lead to the wrong description of 
orbital-ordering with standard DFT methods.®!®] 

Beyond gas-phase molecules, band gaps of various or¬ 
ganic molecular crystals have also been recently success¬ 
fully described with an OT screened-RSH (OT-SRSH) 
functional!®! This was achieved by including a new con¬ 
straint for the asymptotic behavior of the exchange corre¬ 
lation potential, thereby taking into account effects aris¬ 
ing from polarization-induced band renormalization.®*®* 
An open question, however, is whether the OT-SRSH 
functional is capable of accurately predicting not only 
band gaps, but the entire outer valence spectrum of 
molecular crystals. In particular, it is interesting to ex¬ 
amine the OT-SRSH accuracy when the crystal is com¬ 
prised of more complex organic molecules, such as 5Q, 
that are characterized by a mix of localized (<r) and delo¬ 
calized (7r) states as frontier orbitals, where different SIEs 
are expected. It has not yet been investigated whether 
the OT-SRSH approach can accurately deal with self¬ 
interaction problems and at the same time cope with 
polarization effects arising from the crystalline environ¬ 


ment. Capturing both is necessary for an overall good 
description of the electronic structure. For this purpose, 
5Q turns out to be an ideal test candidate. 

In this article, we report a combined experimental and 
theoretical study of the electronic structure of 5Q, which 
answers the above question. We performed angle resolved 
UPS experiments taken on well-ordered films of 5Q in the 
/3-phase, * 49 * 50 * and provide a detailed theoretical assign¬ 
ment of the various peak positions, using the OT-SRSH 
method. To this end, we first investigated the isolated 
5Q molecule, by performing OT-RSH calculations. This 
yielded an IP in excellent agreement with literature data 
from gas phase ultraviolet photoemission spectroscopy 
(UPS). We then calculated the electronic structure of the 
bulk /3-phase, taking into account the correct asymptotic 
behavior in the OT-SRSH approach by computing the 
dielectric constant of the bulk crystal within the random 
phase approximation. To gain a better understanding of 
our results, we also performed many-body perturbation 
theory calculations, within the G 0 W 0 approximatiort®®* 
using various DFT starting points. We obtained excel¬ 
lent agreement with experimental results for both the 
OT-SRSH and G 0 W 0 calculations, the latter based on a 
DFT starting point obtained from a conventional hybrid 
functional. 


II. EXPERIMENT 

5Q films were grown in-situ in ultra-high vacuum 
(UHV) at room temperature on an atomically clean and 
ordered Cu(110) substrate. The Cu surface was prepared 
by repeated cycles of Ar + -ion bombardment and anneal¬ 
ing at 800 K. A 260 A thick 5Q film was deposited in 
situ from a thoroughly degassed evaporator, such that 
the pressure in the system remained at the 10 - 10 mbar 
range during film growth. The nominal growth rate was 
2 A min -1 , as monitored by a quartz microbalance as¬ 
suming a density of 1.47 g cm -3 . 

The electronic structure has been characterized in-situ 
with UPS. Angle-resolved photoemission (ARUPS) ex¬ 
periments were performed using a VG ADES 400 spec¬ 
trometer described elsewhere! 53 The ADES system was 
equipped with a noble gas discharge lamp (unpolarized 
helium I radiation, hv = 21.2 eV) and a movable electron 
energy analyzer, allowing angle resolved ultraviolet pho¬ 
toelectron spectroscopy in the specular plane, with an 
angular resolution of 1° and a total energy resolution of 
150 meV at room temperature. A photon incidence angle 
of a = 60° was used. After the ARUPS measurements, 
the sample was removed from vacuum for geometrical 
structure investigations, ex-situ, using x-ray diffraction 
(XRD) with both 0/20 scans and pole figure analysis 
(Philips X’PERT four circle texture goniometer)!®* 

The XRD data revealed the /3-phase 5Q 

polymorph®^ 

oriented with its (112) plane (blue 
line in Fig. lb) parallel to the substrate, with four 
equivalent domains. As illustrated in Fig. 1, in any one 
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domain of /?-5Q the axes of the molecules are almost 
parallel to the substrate surface (ss 7°). Due to the 
two molecules in the unit cell having their aromatic 
planes at « 70° to each other, and the multiplicity of 
domains, only small angular variations were observed in 
the angle-resolved UPS. Orbital tomography predictions 
of the angular distribution 11 suggested that both n and 
a orbital emissions contribute to the spectra. Thus 
the experimental spectra may be safely related to the 
calculated density of states. 

III. THEORY 

The electronic structure of both the isolated 5Q 
molecule and its /3-crystal structure has been calcu¬ 
lated using two different types of electronic structure 
approaches. The first type is within the framework of 
DFT, where the exchange-correlation potential is ap¬ 
proximated in several different ways: using the gener¬ 
alized gradient approximation (GGA), as parametrized 
by Perdew, Burke and Ernzerhof (PBE)j® The global 
hybrid PBEOf® The short-range hybrid of Heyd, Scuse- 
ria and Ernzerhof (HSE )® 7 And finally the optimally- 
tuned range-separated hybrid (OT-RSH) functional for 
isolated molecule calculations® and the recently pro¬ 
posed optimally-tuned screened range separated hybrid 
(OT-SRSH)® functional for the molecular solid, re¬ 
spectively. The second type of calculation is based 
on many-body perturbation theory, within the GW 
approximation,!® calculated as perturbative ”one-shot” 
G 0 W(P corrections to DFT-based eigenvalue spectra, 
from either PBE or HSE calculations. While the GoWo 
approach has been often applied to molecules and molec¬ 
ular soli ds and is well represented in the literature (see, 
e.g., Refd 7 -SSIIlHss]) ; OT-SRSH based calculations are 
rather new and therefore we provide a concise overview 
of the basic ideas of that approach. For overviews from 
different perspectives, the reader is referred to Refs. 
[251 ES3 EH EH and [68] and specifically for studies of or¬ 
ganic molecular crystals to Ref[45l 

A. Optimally-tuned range-separated hybrids 

The starting point of the range-separated hybrid 
(RSH) c oncept , which is couched within the GKS 
formalism B3ES consists of a partition of the Coulomb in¬ 
teraction asP^ 1 ! 

1 a + f) erf(yr) + 1 — [a + /3 erf( 7 r)] 

r r r 

Here, r is the inter-electron coordinate, erf is the error 
function and a, /3, and 7 are parameters, which in princi¬ 
ple may be freely chosen or determined empirically!® The 
two parts of the split Coulomb operator are treated differ¬ 
ently when computing the exchange interaction. While 


the first term is treated within Hartree-Fock (HF) theory, 
the second one is treated within a standard semi-local (si) 
approximation!® The parameter 7 is the so-called range- 
separation parameter. It controls which of the two terms 
dominates at a given range. It has been repeatedly shown 
that one uniform value for 7 is not sufficient in every 
case, and t hat in fact 7 can be strongly system- and size- 
dependent Therefore, we prefer to determine 

it separately for each system by fulfilling a non-empirical 
condition. This is the point where the optimal-tuning 
comes into play: 7 is chosen such that the difference 
between the energy of the HOMO level and the IP is 
minimized, i.e., we make use of the ionization-potential 
theorem . 19 22 Practically, this is achieved by minimizing 
the following target function?® 

«) = El e W) + IP ^ N + *)] 2 - (2) 

i 

Here, the are ^ ie HOMOs of the (N + i) elec¬ 

tron molecular systems and i is an integer number. The 
IP~ f - a (N + i) are the corresponding ionization potentials, 
which are calculated from the total energy difference be¬ 
tween the N+i electron and the N+i-1 electron system. 
For the moment, a remains as a free parameter and the 
tuning is performed for each choice of a separately, yield¬ 
ing different optimal 7 values that minimize J 2 ( 7 ; a). In¬ 
cluding i = 0 in the sum of Eq. ([ 2 ]) equalizes the HOMO 
of the neutral system with the IP, while for i = 1 the 
IP of the anionic system, i.e., the electron affinity of the 
neutral system is considered, such that its difference from 
the LUMO eigenvalue is minimized!® By extension of the 
sum to negative values of i, states beneath the HOMO 
are also accounted for. This has been shown to assist 
in obtaining optimal tuning if the corresponding orbitals 
are strongly localized!® 

The remaining parameters a and /3, appearing in 
Eq. 0 , control the behavior of the Fock term at its lim¬ 
its. It tends to f when r —> 0 and to when r —> 00 . 
The asymptotic behavior of the xc functional for r —> 00 
has been shown to be crucial for obtaining accurate gaps 
between the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO) 
and introduces a second constraint on the parameters. 
For an isolated molecule, the correct - r asymptotic be¬ 
havior is thus achieved by enforcing a + j3 = lP 3 As 
a consequence, the senri-local contribution in the long 
range is set to zero and a now controls the amount of 
nonlocal Fock exchange in the short range. This can be 
seen in the expression for the exchange-correlation energy 
of the RSH functional, obtain with this constraint!® 

E xc = (1 — a )Esi R x + aEfj R x + E^ R X + E s i tC , (3) 

where si, x and si, c denote semi-local exchange and cor¬ 
relation, respectively, and HF, x denotes non-local Fock 
exchange. When moving from an isolated molecule to an 
organic crystal, the asymptotic behavior of the xc poten¬ 
tial must take into account dielectric screening effects in 
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FIG. 1. (a) Schematic view of the 5Q molecule, (b) Alignment of the 5Q molecules in the /3-crystalline structure, with the 

(112) plane indicated as a blue line. 


the bulk. Thus for r —► oo, the correct limit should be 
— where s is the scalar dielectric constant. This can 

er ’ 

be achieved by choosing a and /3 such that the condition 
a + /3 = j is fulfilled. Note that a gas phase calculation 
may be seen as a special case of this constraint with e 
being equal to one. For the case e ^ 1, the expression 
for the exchange-correlation energy becomes: 

E xc = (1 - a )Esi R x + aE H I F,x + (1 _ ~) E si R x+ 

+ -E^ x + E s i jC . (4) 

Comparing Eqs. ([3]) and Q, clearly the new condition af¬ 
fects only the LR part of the exchange correlation term, 
i.e., the LR parts get screened appropriately as /3 changes 
from 1 — a to \ — a. With the constraints introduced so 
far, there is still no unique choice of a. Different methods 
based on first principles considerations have been sug¬ 
gested to overcome this deficiency. One of the methods 
is based on a further property of the exact KS potential, 
namely the piecewise linearity of the total energy with re¬ 
spect to the (fractional) particle number. For example in 
Refs. 1701 and 1551 a was chosen such that the curvature, 
and therefore the deviation from linearity, of the total 
energy versus particle number curve is minimized. Ac¬ 
cording to Stein et al.f® it may be possible to obtain the 
optimal a directly from minimization of the target func¬ 
tion J 2 (y; a), because deviations from piecewise linearity 


and from the IP theorem are two sides of the same coin. 
For cases where this fails to indicate a unique optimal 
value of a, a more pragmatic approach which uses the 
good agreement of shif ted PBEO results with experiment 
has been suggested!® In that approach, a is determined 
so that the energy difference between the highest occu¬ 
pied delocalized state and the highest occupied localized 
state best correspond to a reference PBEO calculation. 


B. Computational details 

All geometry optimizations have been performed using 
the PBE functional. In order to circumvent issues con¬ 
cerning van-der-Waals interactions, which are poorly de¬ 
scribed in standard GGA and hybrid functionals,® ESI we 
employed the empirical correction scheme of Grimrne 79 
during the geometry optimization of the bulk struc¬ 
ture. Note th at we have taken lattice parameters 
from experiment! 49 ! 50 ! and only optimized the inter¬ 
nal degrees of freedom. The electronic structure of 
the isolated quinacridone molecule were obtained using 
QCHEM version 4.d® with the cc-PVTZ basis set!® 
DFT solid-state calculations of the crystal /3-phase were 
performed using the PARATEC planewave package)® 
modified to include the new SRSH functional.®! 
Within PARATEC, GGA-based Troullier-Martins norm- 
conserving pseudopotentiald® were employed® to repre- 
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sent the core electrons and nuclei. We used a Monkhorst- 
Pack grid of 3 x 3 x 2 /c-point^S. 

Perturbative GoWo results, including the DFT results 
that serve as their starting point, were obtained using the 
VASP packagefSMSSlwith both the PBEPand the HSEP 
functionals. The projector augmented wave (PAW}p! 
approach was employed for treatment of the core elec¬ 
trons, allowing for a relatively low kinetic energy cut-off 
of about 400 eV. The same 3x3x2 k Monkhorst-Pack 
grid employed above was used. 

Go Wo calculations were performed with full-frequency 
integration, using 48 frequency grid points and approx¬ 
imately 4000 unoccupied states to obtain the dielectric 
function and the self-energy. The scalar dielectric con¬ 
stant was determined as | of the trace of the dielectric 
tensor obtained within the random phase approximation, 
including local field effects.^ 


IV. RESULTS AND DISCUSSION 
A. Gas phase of quinacridone 

In this section, we present theoretical results for an 
isolated 5Q molecule, in order to characterize its elec¬ 
tronic structure and frontier orbitals. Importantly, these 
calculations are crucial for obtaining the optimal a and 
7 parameters, used in subsequent bulk calculations. The 
optimal-tuning process, following Eq. |2]), has been de¬ 
termined using both i = 0,1 and i = —1,0,1. Both 
calculations yielded similar results for the optimized val¬ 
ues of 7 , which deviated from each other by no more than 
0.004 Bohr -1 . All results given here are the results of the 

1 = 0,1 tuning. 

Fig. 2a shows the eigenvalue spectra of an isolated 5Q 
molecule calculated with various xc functionals: PBE, 
PBE0, and OT-RSH results with different values for the 
parameter a. For each value of a and the correspond¬ 
ing p=1-a, the range-separation parameter 7 has been 
optimized separately. The PBE calculation results in a 
HOMO (red) value of 5.0 eV and a HOMO-LUMO (blue) 
gap of only 1.8 eV. The former is too low, by more than 

2 eV, compared to the gas phase UPS IP value of 7.23 
eV (shown as a red dotted line in Fig. 2 )PHln the PBE0 
calculation, in which a 25% of non-local Fock exchange is 
included, the HOMO level is 5.8 eV, still showing a large 
underestimation of the experimental reference, and the 
band gap increases to 3.3 eV. 

Next we consider the 5Q HOMO and LUMO values 
obtained with the OT-RSH approach. We have varied a 
from 0.0 to 0.5 in steps of 0.1 and optimized 7 for each 
a value. We observe that the HOMO value is improved 
to ~ 7.15 eV, which slighlty increases for increasing a 
(within ± 0.05 eV), in excellent agreement with the ex¬ 
perimental value of 7.23 . 91 The HOMO-LUMO gap is 
increased to approximately 6.1 eV, and is only slightly 
(± 0.05 eV) affected by the choice of the parameter a. 

Turning to the outer-valence spectra, PBE yields an 


orbital ordering of the lower-lying occupied states that 
is different from that of all other calculations. In par¬ 
ticular, the states shown as orange and yellow lines in 
Fig(2} i.e., the HOMO-1 and HOMO-2 of the PBE calcu¬ 
lation, are located at higher energies relative to all other 
orbitals. The PBE0 results show a change in the eigen¬ 
value spectra, along with energy-level ’’stretching”. In 
our OT-RSH results, the orbital ordering is parameter- 
dependent: all orbitals with 7r-symmetry (which clearly 
exhibit a similar degree of delocalization) show only lit¬ 
tle sensitivity to a. In contrast, the er-orbitals (HOMO-1 
and HOMO-2 of the PBE calculation) are strongly af¬ 
fected and are shifted downwards by more than 1 eV 
when a changes from 0.0 to 0.5. This is not surprising: 
similar observations were made by Refaely-Abramson et 
alP^ and confirmed in additional studies .HHMM] 


In order to provide an explanation for the origin of 
this behavior, orbital plots of the five highest-occupied 
orbitals, as well as the LUMO, are shown in Fig. 2b. 
Note that the color code for the isosurfaces is the same 
one used in the level diagrams of Fig. 2a. Comparing 
the shape of all orbitals in the probed energy range, one 
recognizes the a symmetry and the higher degree of local¬ 
ization of the two orbitals mentioned ab ove. Building on 
experience with other organic molecules j 38 l 40 ‘H S|9ll9 7 the 
reason for the different orbital energies and ordering be¬ 
tween PBE and hybrid calculations is assigned to their 
different self interaction error (SIE). It should also be 
noted that these two cr-orbitals are the main difference 
between 5Q and pentacene, which has no cr-or bitals in 
the energy range of 5 eV below the HOMO j 41 l 42 l 98 l It has 
been previously shown 49 that all outer-valence frontier 
orbitals of pentacene exhibit similar SIE, and therefore 
orbital ordering in pentacene is less sensitive to the choice 
of the xc functional. 


Finally, we also observe that the optimized 7 parame¬ 
ter decreases with increasing amount of short-range Fock 
exchange. This can be rationalized by the range 4 at 
which full Fock exchange sets in, which can be extended 
to larger distances if the a mount of Fock exchange at SR, 
governed by a , is increased! 38 * 39 * 92 lFig. 2 c shows the min¬ 
imal J 2 , obtained for the optimized 7 value for each a, as 
a function of a. The curve shows a distinct minimum of 
J 2 ( 7 ; a) for a val ues be tween 0.2 and 0.3 (note the scale 
bar). It was showrP^Il that there is a rigorous quantita¬ 
tive equality between deviations from piecewise linearity 
and deviations from the IP theorem, represented by J 2 . 
We therefore chose the a value of 0.25, that minimizes J 2 , 
to study the electronic structure of the bulk. This opti¬ 
mal value remains unchanged when including z = — 1 , 0,1 
in the 7 -tuning, or when comparing the energy difference 
between localized and delocalized states with PBEOp^as 
discussed above. 
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FIG. 2. (a) DFT eigenvalue spectra of an isolated 5Q molecule, calculated with PBE, PBEO, and OT-RSH for for various 

values of the SR Fock-exchange parameter a. For each choice of a, the optimal range-separation parameter 7 is also denoted. 
The same color represents the same orbital in the different calculations. The red-dotted line denotes the experimental gas- 
phase UPS result for the IP. 9 h (b) Images of the orbitals in the same color code, with the ordering resulting from the OT-RSH 
calculations with a = 0.25. (c) Deviation from optimal tuning, J 2 ( 7 ; a) as a function of a, using the optimal value for 7 in 
each case. 
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FIG. 3. Band gaps of the 5Q /3-crystalline structure obtained 
from different calculations. The peak maxima of the HOMO 
and LUMO derived bands in the corresponding density of 
states curves are shown as horizontal lines. In each case, the 
middle of the band gap is aligned with that of the computed 
gas-phase HOMO-LUMO gap. 


B. Solid /3-phase of quinacridone 

We now turn our attention to 5Q in the solid /3- 
phase. According to Eq. Q, the scalar dielectric con¬ 
stant £ governs the asymptotic behavior of the xc po¬ 
tential. For organic molecular crystals, the short-range 
interactions are mainly governed by the molecule prop¬ 
erties. We therefore take a = 0.25 and the optimized 
7 = 0.15, as obtained from the above-discussed OT-RSH 
calculation for the isolated molecule. We take £ to be 
3.5, a value obtained from RPA calculations based on 
PBE eigenvalues, because it is already available as a by¬ 
product of our G 0 W 0 results. Note, however, that it 
could easily b e ta ken from computationally inexpensive 
approachesJUllilBJ anc j that this, in fact, is the recom¬ 
mended procedure if a comparison with GW is not per¬ 
formed. 

The results of the various DFT and GW calculations 
for the band-gap are summarized in Fig. 3. Note that 
due to the crystal structure exhibiting two molecules per 
unit cell, each molecular state splits into two bands in 
the crystal. In Fig. 3 we have defined the band gap 
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as the peak-to-peak energy difference derived from the 
computed density of states. Also note that in the bulk 
calculation the absolute energy position of the highest- 
occupied and lowest-unoccupied levels are ill-defined, due 
to the lack of a reference vacuum energy. Therefore, we 
have aligned the center of the band gaps with those of 
the corresponding gas phase calculations. It must also be 
noted that, strictly speaking, a comparison of the com¬ 
puted ’’bulk” IPs and EAs with experimental values de¬ 
termined from thin films is problematic due to surface 
effects not being accounted for in the calculation ! 48 ! 104 ! 

When comparing the values for the fundamental gaps 
we find that, as expected, the PBE gap of about 1.6 
eV is much smaller than those obtained with all other 
approaches. More importantly, it is also approximately 
unchanged compared to the isolated molecule, because 
no effects arising from the polarization of the env iron- 
ment are accounted for with the PBE xc functional .^ 4 ^ 
When computing GoWo corrections to the PBE eigen¬ 
values (GW@PBE), the gap increases to 3 eV. Note that 
polarization effects are inherently taken into account in 
these results as they are contained in the self-energy ex¬ 
pression, leading to a considerably smaller gap than for 
the isolated moleculeWith our OT-SRSH calculation, 
we obtained a band gap of 3.1 eV, which is essentially 
the same as for GW@PBE. Thus the bulk band gap is 
roughly halved, compared to the corresponding calcula¬ 
tion of the HOMO-LUMO gap in the isolated molecule. 

While for the isolated molecule the choice of a barely 
influences the band gap, for the bulk we observe a slightly 
different behavior. When reducing a from 0.25 to 0.0, but 
keeping e = 3.5 unchanged, we observe a 0.3 eV reduc¬ 
tion of the band gap. This finding can be explained by 
considering Eq. ©■ When reducing a , the optimized 7 
increases. In other words, when decreasing the amount 
of SR Fock exchange, the amount of LR Fock exchange 
increases to maintain the amount of overall non-local ex¬ 
change. As a consequence, the spatial LR region in which 
the effective dielectric screening acts extends, thereby en¬ 
hancing polarization effects and leading to a smaller band 
gap. 

Now we compare the PBE results with those based on 
a HSE calculation. HSE yields a band gap of 2.3 eV, 
which lies between the pure PBE and the corresponding 
OT-SRSH result. In order to provide an explanation, 
recall the properties of the HSE functional. As already 
mentioned, it is a SR hybrid functional using non-local 
exchange only in the SR and pure semi-local exchange 
in the LR. The amount of Fock exchange is given by 
a = 0.25 and a universal 7 value of 0.11 is used. Thus, 
the result is improved over PBE as some non-local ex¬ 
change is introduced. However, because the xc potential 
decays exponentially, i.e., e = 00 in Eq. 0 . the asymp¬ 
totic behavior is incorrect and a smalle r ban d gap than in 
the OT-SRSH calculation is obtained!^ 1 ®! Compared to 
PBE, the increased HSE band gap results in the fact that 
a subsequent RPA calculation yields a decreased scalar 
dielectric constant of e = 2.8. When using this value 


in the OT-SRSH calculation, we obtain a band gap of 
3.3 eV, which is larger compared to the OT-SRSH band 
gap achieved using the PBE-based RPA e due to the re¬ 
duced screening with the smaller HSE-based e. Again the 
band gap compares well with a Go Wo computation with 
an HSE starting point (GW@HSE), which yields a band 
gap of 3.5 eV. The SRSH gaps are then indeed consistent 
with the GW calculation, given a similar scalar dielectric 
constant. 



FIG. 4. Comparison of experimental and theoretical pho¬ 
toemission spectra for the /3-phase crystalline structure of 
5Q. (a) Experimental UPS data in normal- and 50° emis¬ 
sion (black solid and dashed lines, respectively), (b) Com¬ 
puted results of PBE and GW@PBE (blue dashed and solid 
lines,respectively), (c) Computed results of OT-SRSH ap¬ 
proach with a = 0.0 and a = 0.2 (red dashed and solid lines, 
respectively), (d) Computed results of HSE and GW@HSE 
(orange dashed and solid line, respectively). For all spectra 
shown the energy axis has been aligned with the maximum 
of the highest occupied peak, which has been set to zero. A 
Gaussian broadening of 0.2 eV was used in all computed spec¬ 
tra. 

































Next, we investigate the valence band electronic den¬ 
sity of states of the crystal, calculated at various levels 
of DFT and G 0 W 0 , in more detail. These are compared 
with experimental ARUPS results. The results are sum¬ 
marized in Fig. 4. Note that both experimental and cal¬ 
culated spectra have the energy axis shifted such that 
the highest occupied peak maximum has been aligned to 
zero. A Gaussian broadening of 0.2 eV was used in all 
computed spectra. 

In the experimental spectra, there are four peaks in 
the shown energy range, which are highlighted by ver¬ 
tical dashed black lines. The first peak coincides with 
0 eV, by virtue of the alignment procedure. The other 
peaks are located at —1.1, —2.0 and —3.4 eV. For fur¬ 
ther interpretation of the experimental findings and for 
gaining more insight into the origin of individual peaks, 
we compare the peak positions and spectral shape of the 
experimental data to the theoretical results. In Fig. (4 Jd, 
we display the PBE result (blue, dashed line) and the 
corresponding GW@PBE calculation (blue, solid line). 
At first sight, agreement between the PBE result and 
experiment appears reasonable because the peak posi¬ 
tions seemingly agree quite well. However, when taking 
into account the peak shape as well, this agreement turns 
out to be coincidental. As in the case of the isolated 
molecule, PBE provides the wrong orbital energies and 
ordering, particularly for the localized states of a symme¬ 
try. The GW@PBE results certainly improve the band 
gap, as shown in Fig. [3j and also somewhat stretch the 
valence band spectrum. However, judging by the com¬ 
parison with the experimental data, the GW@PBE result 
does not seem to accurately describe the investigated sys¬ 
tem. The main deviation of the GW@PBE curve is its 
peak at —1.4 eV, which does not show up in the experi¬ 
mental data at all. 

In order to clarify this deficiency of the GW@PBE re¬ 
sult, we compute the valence band DOS resulting from 
the OT-SRSH approach using a = 0.25, shown as a red 
solid line in Fig. 4c. From the previous findings for the 
isolated 5Q molecule, as well as similar molecules,® we 
know that the amount of SR Fock exchange mainly af¬ 
fects states with a distinct degree of localization. Such 
behavior is also expected for the bulk. Therefore, Fig. & 
shows in addition the DOS obtained from an OT-SRSH 
calculation with a = 0.0 as a red dash-dotted line. Indeed 
the two mentioned spectra in Fig. [4J: are dramatically dif¬ 
ferent, although only little influence of a on the size of the 
band gap was observed for the molecule (Fig. [2]) and the 
molecular crystal, as discussed above. Interestingly, the 
a = 0 curve resembles the GW@PBE result, including a 
peak at about —1.4 eV, while the a = 0.25 spectrum has 
no peak at that energy, and the corresponding states are 
shifted to lower energies. This enhances the peak at —2 
eV and leads to a rather impressive agreement with the 
experimental data. Thus, the origin of the incorrect peak 
at —1.4 eV in the GW@PBE and the OT-SRSH with 
a = 0 is related to a remaining SIE of the strongly local¬ 
ized cr states (depicted in yellow and orange in Fig. [2|, 


resulting with calculated energy levels that are too high 
for these states and changing the overall spectral shape. 
This is confirmed by plots of the orbital density associ¬ 
ated with these states. In Fig. [5] a density plot of the 
HOMO-3 orbital of the isolated molecule obtained from 
OT-RSH with a = 0.25, as well as the partial charge- 
density of the corresponding orbital in the solid-state, is 
shownP® The extended bulk-state can be clearly associ¬ 
ated with the respective orbital of the isolated molecule 
and the same is true for all other outer valence states 
(not shown for brevity). 



FIG. 5. (a) Density plot of the HOMO-3 orbital obtained 

from the OT-RSH calculation of the isolated molecule, (b) 
Partial charge density decomposed on the band derived from 
the orbital shown in (a). 

Finally, we compare these results with those based on 
a HSE calculation. The spectra of the HSE and the cor¬ 
responding GW@HSE calculation are shown in Fig. [4jd as 
orange dashed and solid lines, respectively. Moreover, the 
OT-SRSH spectra with a = 0.25 and using e = 2.8 from 
the HSE-based RPA calculation is displayed in Fig. 4c 
as a red dashed line. Clearly the two SRSH calculations 
with a = 0.25 coincide, up to a very small deviation. 
This observation shows that the screening introduced in 
Eq. Q affects all occupied states similarly and shifts 
the whole spectrum rigidly, thereby changing the com¬ 
puted band gap appropriately. Namely, while the e value 
greatly affects the gap renormalization (as it is the mea¬ 
sure of electrostatic polarization), the shifted occupied 
spectra is almost entirely dependent on the values of a 
and 7 , and is practically the same for the two examined 
values of e. The line shape of the outer-valence band 
spectrum of the HSE calculation is in good agreement 
with experiment and OT-SRSH calculations (a = 0.25). 
The GoWo calculation with the HSE starting point shift 
the HSE spectrum rigidly and we find an almost per- 
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feet agreement with the experiment. Furthermore, it be¬ 
comes obvious that the full-frequency GW@PBE calcu¬ 
lation suffers from the inappropriate starting point, as 
already reported in a number o f stu dies of the organic 
and metal-organic molecule.— 107 1!2SJ 

Having found the theoretical methods which yield an 
accurate description of the electronic structure of the or¬ 
ganic molecular crystal, we are able to assign specific 
molecular states to the experimentally observed peaks 
(a band-decomposed charge-density plot of one of these 
states was given in Fig. §>)■ In Fig. [6] we follow such 
plots to assign molecular orbitals to the appropriate peak 
in the OT-SRSH functional calculation with a = 0.25. 
This optimal fraction of SR Fock exchange allows a si¬ 
multaneous prediction of both cr-type localized orbitals 
(yellow and orange in Fig. [6| and 7r-type delocalized or¬ 
bitals (other orbitals in Fig. |6j) , as discussed above. By 
that, it allows the assignment of theoretical orbitals to 
peaks of the experimental UPS data, as shown in Fig. [6] 



E ‘ E HOMO < eV ) 


FIG. 6 . Measured (normal and 50°) photoemission spectra 
for the /9-crystalline structure of 5Q, including density plots 
of the contributing molecular states, as well as the energy 
ordering shown as vertival lines, obtained using the OT-SRSH 
functional with a = 0.25, 7 = 0.15 and e = 2.8. 


V. CONCLUSION 

In conclusion, we have studied the electronic structure 
of the organic molecule quinacridone in the gas phase and 
in the crystalline /3-pliase. For the gas phase, we find that 
the semilocal PBE and the hybrid PBEO functional both 
underestimate the ionization potential and the funda¬ 
mental band gap. When employing an optimally-tuned 
range separated hybrid functional, we obtain an excellent 
agreement with the experimental UPS value for the ion¬ 
ization potential and a larger fundamental gap, which is 
insensitive to the fraction a of short-range Fock exchange. 


The latter parameter, on the other hand, is shown to be 
crucial for attaining the correct relative orbital ordering 
of delocalized versus localized orbitals. 

For the crystalline bulk /3-phase of quinacridone, 
we have conducted angle-resolved photoemission experi¬ 
ments, which were used to benchmark our calculations. 
We demonstrate that by using the optimally-tuned value 
for the fraction of short-range exchange (a = 0.25), the 
correct orbital ordering was obtained within the occupied 
manifold of the states. We further showed that in order to 
take into account the level renormalization due to elec¬ 
tronic polarization in the crystalline phase, the appro¬ 
priate asymptotic behavior of the exchange-correlation 
functional is essential. This has been achieved by employ¬ 
ing the optimally-tuned screened range-separated hy¬ 
brid (OT-SRSH) approach, in which the screening is ac¬ 
counted for by an effective dielectric constant, e, which 
we computed from the trace of the macroscopic dielec¬ 
tric tensor obtained within the random phase approxi¬ 
mation. Thus, the description of the unoccupied mani¬ 
fold of the states, in particular the fundamental gap, is 
also greatly improved, showing a band gap renormaliza¬ 
tion from the gas phase to the bulk based on physical 
grounds. For comparison, we computed the GqWq cor¬ 
rected electronic structure of the bulk using both PBE- 
GGA and the short-range hybrid functional HSE. These 
results emphasize the importance of the starting point 
in this perturbative approach, where the GoWo@HSE es¬ 
sentially agree with the OT-SRSH approach. 

In summary, our work shows that the reliability of the 
optimal-tuning approach for molecular systems can be 
extended to the valence spectrum of molecular solid sys¬ 
tems and that results at a level of accuracy comparable 
with GW calculations can be achieved. We emphasize 
that, based on physically motivated choices for the pa¬ 
rameters, the OT-SRSH approach allows for an accurate 
description of the band gap and at the same time of the 
relative orbital energies of the outer valence spectrum 
without any empiricism. Therefore, it may serve as a 
computationally inexpensive and reliable tool. 
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